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Abstract 

An exact formula for the minimal coordination numbers of the parallel packed bundle of rods is 
presented based on an optimal thickening scenario. Hexagonal and square lattices are considered. 

PACS numbers: 87.15.-v, 87.16.Ka, 61.46.Fg 



I. INTRODUCTION 



Hexagona 



DNA 
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packing of slender objects (filaments or rods) like semiflexib 

26| . collagen 



15| . filamentous actin 



12| or nanotubes [5|, llSl . 



e p olymers, e.g. 



25| provides the 



most compact way of filling a domain in space. It maximizes the number of contacts between 
the rods. It has long been known that DNA can condensate to form bundles, in particular, 
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17( 1 . Similar 



closed toroids when polyvalent cations are added to the solution [sl. Isl. [lO. 
shapes are observed for other filaments stabilized by different types of adhesion: by hydrogen 
bonding, by ionic interactions or by van der Waals forces ly, Q]- ^ perfect hexagonally 
packed bundle must be twistless to provide continuous parallel arrangement of filaments 



though it can writhe in space 



22|. On the other hand, the interaction peculiarities of 



individual filaments in a parallel bundle may favour them to fill in a non-hexagonal (e.g., 
square) 2D lattice in cross section. 

The main purpose of this paper is to present a compact exact formula for a number of 
unoccupied bond sites on the surface of a tightly packed parallel bundle. This quantity 



13, 19, 20, 21 and it has 



commonly appea. in analys. of ^tabjHtvof DNA to.oid. fl fl 

been computed either approximately [2, l20| or by a complicated procedure [19|, |23| or for 
particular numbers of strands in the bundle 0, [l^ . It is reasonable to believe that the exact 
analytical expression will prove to be handy in the future analytical research. 

Consider a bundle of parallel packed rods. In the Euclidean space, each rod may be in 
continuous contact with at most 6 other rods [21]. The total number of binding sites is 6A^ 
for hexagonal lattice and AN for square one. Each actual bond involves two sites. Therefore, 
the number of the binding sites that are exposed to the environment is even and one half of 
it is called the coordination number 201 ]. 



In a poor solvent, there is a relative energetic advantage for a filament to have a polymer- 
polymer contact versus being exposed to solvent. Therefore, we are interested in the tightest 
bundle when the coordination number is minimal. The minimal coordination number enters 
a free energy expression that is used in studies of DNA condensation j^, 4 
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20 



23]. 



Its behaviour as a function of number of filaments affects stability of various conformations. 

Since the rods are parallel, all orthogonal cross sections of the bundle are congruent and 
we can treat the problem reduced to two dimensions, i.e. we consider the set of discs that 
appear in an orthogonal cross section of the bundle. The minimal coordination number 



2 



equals 



for the hexagonal lattice and 



OlN 



^Jl2N - 3 



N 



2VN 



(1) 



(2) 



for the square lattice. In Eqs. ([I]),(l2l), [■] stands for the ceiling function: for a; G M, [x] is 
defined as the smallest integer greater than or equal to x. In particular, ai = 3, a2 = 5, 
as = 6, and so on and Pi = 2, P2 = 3, = 4, and so on. 

Equation ([1]) gives exactly the same values (up to factor of 2) as a rather overcomplicated 
algorithm that was suggested in Ref. 3], Eq. (3), and was used recently in Ref. [2^. An 
approximate expression was presented earlier for the average number of interactions which 
equals 6 — 2Q;Ar/A^ in our terms (see Eq. (10) in Ref. j^). That approximate formula may be 
used to compute the coordination number only for small number of members, up to < 30. 
Actually, Eq. (10) describes a nonmonotonic function which contradicts the physical sense. 
Note also that the graph in Fig. 3 in the cited article which is claimed to show the function 
from Eq. (10), in reality, represents a different non-specified function. 

For a particular case of the bundles with a complete hexagonal cross section, one can find 
an exact expression for the average number of interactions in Ref. {4] (Eq. (4)) which agrees 
with Eq. a N = 3n'^ + 3n + l,n = 1,2, . . .. It was suggested to use the formula of the 
„u,„be.. of bo„d. computed fo. complete hcxago,. to approximate the opt.mal coufigu.afo,. 



for artbitrary number of members 13|]. The limiting coordination number = 2v3A^ 
proposed in Ref. [2^ approximates Eq. ([T]) for a hexagonal bundle of A^ filaments when 
N ^ 00. No proof was given on optimality of conformations considered in the above 
mentioned works to compute the minimal coordination number. 



II. PROOF OF EQ. 

1. We show that Eq. ([1]) is valid for complete shells, i.e. for A^ = M{n) = 3n'^ — 3n + 
l,n = 1,2, .. . (Fig. [1]). Clearly, ai = 3 for n = 1. For n = 2,3, . . ., there are 6 corner 
cells, each with 3 open bonds and Q{n — 2) side cells, each with 2 open bonds. In total, 
6 X 3 + (6n - 12) X 2 = 2(9 + 6n - 12) = 12n - 6 = 2ajv/(„,) Q. 

2. It is easy to check the correctness of Eq. ([T]) for A^ = 1,...,M(2) = 7 by direct 
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FIG. 1: Complete shells of generation n = 1,2,3. Here and in all other figures, bonds are shown 
as black sticks, free potential bonds are dashed. 

computation (Fig. [2]). 





FIG. 2: Optimal configurations for = M(l) = 1, . . . , M(2) = 7. 

Note that, as shown in Fig. [3], the configuration with minimal coordination number is 
obhged to be neither unique nor convex. (Convexity may be defined as absence of cells with 
only one free bond.) 

3. We show that Eq. ([T]) is valid for = M(n), M{n) + 1, . . . , M{n+ 1) for any generation 
number n > 2 for a chosen algorithm of adding cells. 

Consider what happens when we add more cells to the complete shell of generation n 
(n = 3 in Fig. Hj). The first added is bound to two other cells independently of its place 
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FIG. 3: All three configurations have the minimal coordination number oq = 9. 

(A^ = 20). So, it closes 2 existing open bonds and adds 4 new which means that the 
coordination number increases by 1: Qm+i = c^m + |(4 — 2) = + 1- To minimize the 
number of exposed sides, the next cell has to be added adjacent to the first one so that 
it contacts with three cells (A^ = 21). Since the number of new open bonds is 3 as well, 

We can keep adding new cells until the entire side of the hexagon is filled up. As it 
happens, the coordination number increases by 1, because the (M(n) + n)-th cell can only 
be bound with two others (A^ = 22). The cell has to be added to one of the remaining 5 sides 
of the initial complete shell. One of the neighbouring sides has to be chosen to make the next 
increase of the coordination number as late as possible (A^ = 23). As the neighbouring side 
is being filled, the coordination number does not change until the entire row is completed 
(A^ = 24). The cycle repeats six times until all six sides of the {n + l)-th layer are filled up 
(only filling two more sides is shown in Fig. H]). Each time as the (M(n) +m)-th cell takes its 
place [i = 1, 2, 3, 4, 5), the coordination number adds 1. There is no need for rearrangement 
of cells as more cells being added. Summing up, we can say that the coordination number 
increases each time when a cell is added to a convex configuration. 

Let us now check whether Eq. ([1]) follows precisely this pattern. As seen from the equation, 
On = ttAf-i + A, for integers A^ > 1 and A > 0, if and only if there exists an integer K such 
that 12(A^ -l)-3< K"^ < 12N - 3. Now check what are A^'s for which A > 1. The latter 
would imply ^/l2iV^^ — a/12(A^ — 1) — 3 > 1 which only possible for A^ = 2, 3. Since we 
are now interested in the case A^ > M(2), we are safe to consider only A = 1. The changes 
in may be indexed with K. 

There are 6 changes in a^v between the shells of two consecutive generations, i.e. in the 
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N=20 N=21 N=22 N=23 




N=24 N=25 N=26 




N=27 N=28 N=29 N=30 




FIG. 4: Adding more cells to the complete shell of generation n = 3. A choice of optimal sequence 
for N = 20,... , 30. 

range of numbers between M{n) and M{n + 1). Showing first that = + 1, let 

us set K = aM{n) = 6n - 3. Then l2M{n) - 3 < (Gn - 3)^ < 12(M(n) + 1) - 3. Indeed, 
12(3^2 _ 3n + 1) - 3 = (6n - 3)^ < 12(3n2 _ 3^ + 2) - 3. 

As explained above, the remaining 5 increases oia^ should take place for = M(n) +in, 
2 = 1, 2, 3, 4, 5. To prove this, take K = aM{n) + h then we can write 

12(M(n) + m - 1) - 3 < (ttMH + i? < 12(M(n) + in) - 3, 
12(3^2 - 3n + 1 + m - 1) - 3 < (6n - 3 + z)^ < 12(3n^ - 3n + 1 + m) - 3, 
36^2 + 12(i - 3)n - 3 < 36n2 + 12{i - 3)n + (i - 3)^ < 36n^ + 12(i - 3)n + 9, 

which is obviously satisfied. Therefore, aM(n)+in = C(M(n)+in-i + 1- 

Thus we see that, according to our formula, jumps by 1 when the first cell is added to 
the complete shell and when a portion of n cells is added to the complete shell of generation 
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n unless the n + 1 layer is completed. By way of induction we conclude that Eq. ([T]) is valid 
for all > 1. 

4. The remaining thing to do is to show that the above scenario is optimal in the sense 
that it provides the minimal number of exposed sites for every number of cells. Clearly, an 
optimal configuration must be simply connected. Let us take an arbitrary simply conected 
cluster of cells. There exists a simple closed polygonal line that connects the centres of 
all p cells with exposed bond sites (we call it a perimeter). Let us orient the perimeter so 
that the interior of the cluster (A^ — p cells) be on the left. 

Each perimeter cell can belong to only one of the following 4 groups. The first group 
consists of t cells where the perimeter line makes no turn. Such a cell has two free bonds 
on the right and at most two bonds with the interior part of the cluster. The second group 
includes u cells where the perimeter turns to the left through vr/S. The members of the 
second group have three exposed bond sites on the right and at most one bond with the 
interior. The third group is made up with v cells with the 27r/3 turns of the perimeter and 
with 4 free bonds on the right. The cells of this group are bonded only with two perimeter 
neighbours. The last fourth group is formed with w cells in which the perimeter line turns 
to the right through 7r/3. Such perimeter cells cannot appear in a convex cluster. They 
have one free bond on the right and up to three bonds to the interior. Thus, we have 

p = t + u + v + w. (3) 

When traversing the complete path along the perimeter, we make exactly one turn (the 
winding number of the perimeter is one), hence we can write 

u + 2v -w = (4) 

By definition, the coordination number of the A^-cluster equals half the number of 
exposed bond sites on the right of the perimeter. Summing up these exposed sites gives 

2aN = '2,t + 2,u + 4:V + w. (5) 

Eliminating the quantities t,u,v,w from Eq. ([5]) with the help of Eqs. ([3]), (jl]) leads to 

"AT = P + 3. (6) 

The last equation simply shows that the coordination number is essentially the length of the 
perimeter, i.e. the number of the exposed cells. 
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All the bonds with the interior contribute to the coordination number of a cluster that 
is obtained after peeling off the A^-cluster. Then, the coordination number of this (A^ — p)- 
cluster can be estimated as 

2aN-p <2t + u + 3w. (7) 
Again, with the help of Eqs. ([3]), (jlj) we obtain from Eq. (I7j) 

on-p < P - 3. 

Combining the last inequality with Eq. (jH]) allows us to relate the coordination numbers of 
the full and the stripped clusters 

aN > cuN-p + 6. (8) 

Let us now assume that the coordination number of the peeled-off (A^—p) -cluster satisfies 
Eq. ([T]) for the minimal coordination number. Then we can write 

aN-p > ^/l2{N-p)-3 = ^/12{N - a^) - 33 

(here we used Eq. ([6]) to eliminate p in the right-hand side). Substitution of a^^p from 
the last inequality into Eq. and some algebraic transformations eventually lead to the 
required estimate for the A^-cluster 

On > Vl2N -3. 

Thus we see that it is impossible to lower the coordination number below the value that 
was computed for the sequence of clusters we considered above in the former parts of this 
proof. We conclude that indeed such a scenario of growth guarantees the most compact 
configuration for any number of cells with the minimal coordination number from Eq. ([1]). 

Proof of Eq. ([2]) is similar. 

III. DISCUSSION 

1. The minimal coordination number formula may be used to compute the number of 
actual bonds U{N) = 3N — \\^12N — 3 1 (for hexagonal lattice). For complete hexagonal 
shells with A^ = M{n),n = 1,2,... (Fig. [1]), the ceiling function can be omitted. As 
suggested by Ishimoto and Kikuchi[l3j], the expression V^(A^) = 3N — \^12N — 3 can be used 
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as an approximation for arbitrary A^. They present a comparison of V^(A^) for real with 
what they call "the exactly counted discrete function" (though they plot two continuous 
functions). It is not clear from the paper how "the exactly counted discrete function" was 
computed. One can see from the graph (Fig. 2 in [13]) that the latter function differs from 
U{N) (e.g., notice too big difference between the curves for = 30). Anyway, our theory 
provides a solid proof that V^(A^) can be used as an approximate function for large A^ since 
it cannot differ from the exact optimal function by more than 1. 

2. Equation ([1]) may be "inverted" to give an equation for A^^, the maximal number of 
strands in the bundle for given (integer) coordination number a: Na = [(a^ + 3)/12j , a > 5, 
where [-J stands for the floor function: for x G M, \_x\ is defined as the greatest integer less 
than or equal to x. For a = 9 and a > 11, the A^^-bundles are considered as particularly 
stable in a sense that they have the same surface energy as the {Na — l)-bundles: addition 
of the Na-th filament does not affect the coordination number. However, one more filament 
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20|. 



increases the surface energy. Such numbers Na are called magic and supermagic 
For the square lattice the "inverted" Eq. (El) reads Np = [(3'^ /A\, (3 >2. 

3. A word of warning should be placed here concerning the use of the mininal coordination 
number to estimate the surface energy of the bundle. The coordination number, we consider 
in this paper, is actually defined for a thin cross-sectional slice of a bundle. To compute the 
surface energy of the entire bundle, one need to integrate over its length. Since in general 
the lengths of the particular strands may differ, the total number of the exposed binding 
sites does not necessarily amount to simply a doubled product of the coordination number 
and the length of the bundle axis. To get the correct result one should account for the actual 
shapes of individual filaments exposed to the environment and these shapes could be rather 
complex. This subject is beyond the scope of this paper. 

It is only worth to mention here that, for closed parallel bundles the lengths of the 
components are quantized and this fact simplifies the computation of the surface energy 
reducing it essentially to the coordination number multiplied by the length. However, the 
life is not so simple even in this case because there are intrinsic geometric constraints on the 
arrangements of the rods into the bundle which are not always consistent with the minimality 
of the exposed surface [22j]. It may also be the case that the sequence of optimally packed 
bundles is kinetically unfeasible. For example, when the filament coils into a toroidal bundle 
it is unlikely that it passes through the hole. 
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4. Let us here focus on the hexagonal lattice first. The complete shells are indexed with 
the generation number n and the number of cells is = 3n^ — 3n + l, n = 1,2, . . .. Then, for 
any A^-bundle we can define a real quantity Uai^N) G M, z/„ > 1 such that = 3z/^ — Si/q, + 1. 

The problem of finding a configuration that minimizes the coordination number may 
be considered as an isoperimetric problem where N plays part of the "area content" and 
the coordination number is essentially the perimeter length (see Eq. ([6])). In the classic 
continuum analogue, the solution of the 2D isoperimetric problem is a circle with area 
A = vrr^ and circumference 5* = 27rr so that ^ = S with the radius r being the so- 
called "harmonic parameter" [9(. In the discrete problem, the last equation transforms into 
= [6z/q ^ 3] = |'Vl2A^ — 3 1 = On- Thus, for complete shells = n serves as the 
harmonic parameter and the above equation may be viewed as a discrete analogue of the 
known expression for the perimeter of the regular polygon as a derivative of its area relative 
to the radius of the inscribed circle 31. Note that the latter relationship follows from the 

n 

Steiner-Minkowski theorem Thus, parallels can be seen between the sequence of the 
solutions to the discrete isoperimetric problem and the painting (or collaring) of polygons. 
For the square lattice, we have correspondingly vpiN) = \/N and = /3n- 

5. The proven extremality property of layered packing may look a little counterintuitive 
for thick bundles with a large number of rods. Indeed, one may expect that the optimal 
shape of the cross section would approach a circle as much as possible. As we can prove 
now, this is not the case [2^. Moreover, the relative difference of the minimal shape and the 
circle becomes more noticeable as grows. To see this let us take a complete hexagon with 
M{n) cells. Its perimeter is a regular hexagon and it may be inscribed into a circle of radius 
R = {n — l)d, where d is the diameter of an individual rod. The maximal distance between 
the circle and the inscribed hexagon is proportional to the radius and equals (1 — ^/3/2)R. 
Each of six circular segments appeared between the hexagon and the circle will make a room 
for, say, a couple of additional layers if n > 16 (M(16) = 721). This suggests that these 
gaps could be filled in with two or more layers as more rods join the bundle. On the other 
hand, the optimal scenario of growth prescribes that a layer will be filled completely until 
the next regular hexagon with M{n + 1) is formed. Of course, as we have seen, the optimal 
configurations may not be unique (Fig. [3]) and some of them could be more "circular" than 
others, but diversions to them in the process of growth would imply rearrangements of rods 
which could be prohibited. 
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